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Status of the Standard Model Prediction of the Muon g-2 
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a Dipartimcnto di Fisica, Universita di Padova and INFN, Via Marzolo 8, 35131 Padova, Italy 

The current status of the Standard Model prediction for the anomalous magnetic moment of the muon is briefly 
reviewed and compared with the present experimental value. 
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1. Introduction 

Schwinger's 1948 calculation (]] of the leading 
QED contribution to the anomalous magnetic mo- 
ment of the muon a M = (g^ — 2)/2, equal to the 
one of the electron, was one of the very first re- 
sults of this theory, and one of its early confirma- 
tions. During the last few years, in a sequence of 
increasingly precise measurements, the E821 Col- 
laboration at Brookhaven has determined a M with 
a fabulous relative precision of 0.5 parts per mil- 
lion (ppm) |2I3I4| , serving as an invaluable tool to 
test all sectors of the Standard Model (SM) and 
to scrutinize viable alternatives to this theory 5 . 
This note provides a brief summary of the present 
status of the three contributions into which the 
SM prediction a^ M is usually split - QED, elec- 
troweak and hadronic - and a comparison with 
the current experimental value. 

2. QED and Electroweak Contributions 

The QED contribution to arises from the 
subset of SM diagrams containing only leptons 
(e, r) and photons. The lowest-order contri- 
bution is a« ED (l loop) = a/ (2k) [J. Also the 
two- and three-loop QED terms are known ana- 
lytically - see [2] for an update and a review of 
these contributions. The four-loop term has been 
evaluated numerically, a formidable task first ac- 
complished by Kinoshita and his collaborators in 
the early 1980s [7j. The latest analysis appeared 
in [Hj. Note that this four- loop contribution is 
about six times larger than the present experi- 
mental uncertainty of a M ! The evaluation of the 



five-loop QED contribution is in progress [5j ■ 

Adding up these terms, using the latest CO- 
DATA JU| recommended value for the fine- 
structure constant a" 1 = 137.035 99911 (46), 
known to 3.3 ppb, one obtains [H] a® ED = 
116 584 718.8 (0.3) (0.4) x lCT 11 . The first error is 
mainly due to the uncertainty of the 0(a 5 ) term, 
while the second one is caused by the 3.3 ppb 
uncertainty of the fine-structure constant. 

The electroweak (EW) contribution to a M is 
suppressed by a factor (m^/M^) 2 with respect to 
the QED effects. The one-loop part was computed 
in 1972 by several authors [H]: af w (l loop) = 



5G„m„ 



24\/2ir 2 

where G M = 



l + ±(l-4sin 2 



+ 0{ml/Ml w , H ) , 

1.16637(1) x 10~ 5 GcV~ 2 . Em- 
ploying the on-shell definition sin 2 W = 1 — 
Ml/ 'M 2 [12, where M z = 91.1875(21) GeV and 
M w is the SM prediction of the W mass (which 
can be derived, for example, from the simple 
formulae of jE] leading to M w = 80.383 GeV 
for the Higgs mass M H = 150 GeV), we obtain 
a EH '(l loop) = 194.8 x 10~ n . 

The two-loop EW contribution to a M is not neg- 
ligible because of large factors of ln(M z _ w /m/), 
where mf is a fermion mass scale much smaller 
than M w [HI- It was computed in 1995 [T5] . 
The proper treatment of the contribution of the 
light quarks was addressed in |16I17| . These re- 
finements significantly improved the reliability of 
the fermionic part (that containing closed fermion 
loops) of af^ (two loop) leading, 



150 GeV, to al w = 154(1)(2) x 10" 
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for M H = 

' n P2|. The 

first error is due to hadronic loop uncertainties, 
while the second one corresponds to an allowed 
range of M H £ [114, 250] GeV, to the current top 
mass uncertainty, and to unknown three-loop ef- 
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fects. The leading-logarithm three-loop contribu- 
tion to a^ w is extremely small |17I18| . The results 
of JH| f° r the two- loop bosonic part of a^ w , per- 
formed without the large M H approximation pre- 
viously employed, agree with the previous evalu- 
ation ^S] in the large Higgs mass limit. Work is 
also in progress for an independent recalculation 
based on the numerical methods of |2()j . 

3. The Hadronic Contribution 

The evaluation of the hadronic leading-order 
contribution a" LO , due to the hadronic vacuum 
polarization correction to the one-loop diagram, 
involves long-distance QCD for which perturba- 
tion theory cannot be employed. However, using 
analyticity and unitarity, it was shown long ago 
that this term can be computed from hadronic 
e + e~ annihilation data via the dispersion integral 
< L ° = ^fZl dsK ( s )° m ( s ) EH, where a(°)( S ) 
is the total cross section for e + e" annihilation 
into any hadronic state, with extraneous QED 
corrections subtracted off. The kernel function 
K(s) decreases monotonically for increasing s. 

A prominent role among all e + e~ annihilation 
measurements is played by the precise data col- 
lected in 1994-95 by the CMD-2 detector at the 
VEPP-2M collider in Novosibirsk for the e + e~ — * 
7r + 7r~ cross section at values of y/s between 0.61 
and 0.96 GeV [221 • The quoted systematic error 
of these data is 0.6%, dominated by the uncer- 
tainties in the radiative corrections (0.4%). Re- 
cently [231 the CMD-2 Collaboration released its 
1996-98 measurements for the same cross section 
in the full energy range y/s S [0.37,1.39] GeV. 
The part of these data for y/s <E [0.61, 0.96] GeV 
(quoted systematic error 0.8%) agrees with their 
earlier result published in [22] . Also the SND Col- 
laboration (at the VEPP-2M collider as well) re- 
cently presented its analysis of the e + e~ — * 7r + 7r~ 
process for y/s between 0.39 and 0.98 GeV, with 
a systematic uncertainty of 1.3% (3.2%) for y/s 
larger (smaller) than 0.42 GeV [22- A hint of dis- 
crepancy, at the level of the combined systematic 
error, occurs between the CMD-2 and SND mea- 
surements (the contribution to a" LO of the SND 
data is a bit higher than the corresponding one 
from CMD-2) [23]. Further significant progress is 



expected from the new e + e collider VEPP-2000 
under construction in Novosibirsk |23| . 

In 2004 the KLOE experiment at the DA$NE 
collider in Frascati presented a precise measure- 
ment of er(e + e~ — > 7r + 7r~) via the initial-state 
radiation (ISR) method at the (f> resonance |25| . 
This cross section was extracted for y/s between 
0.59 and 0.97 GeV with a systematic error of 1.3% 
and a negligible statistical one. There are some 
discrepancies between the KLOE and CMD-2 re- 
sults (KLOE's data lying higher than the CMD- 
2 fit below the p peak, and lower on the peak 
and above it), although their integrated contribu- 
tions to a" LO are similar (2H|. The data of KLOE 
and SND disagree above the p peak, where the 
latter are significantly higher. The study of the 
e + e~ — * 7r + 7r~ process via the ISR method is also 
in progress at BABAR [201 and BELLE [23 . On 
the theoretical side, analyticity, unitarity and chi- 
ral symmetry provide strong constraints for the 
pion form factor in the low-energy region |28| . 

The evaluations of the dispersive integral based 
on the analysis [221 01 the 1994-95 CMD-2 data are 
in good agreement: 1 

EH < i0 = 6934 (53) CTp (35) ra(i x 10" 11 , 
[32] al LO = 6948 (86) x 1(T U , 
[33| a H ^° = 6924(59^(24)^x10-", 
EU a% LO = 6944(48) eKp (10) ra(i x 10- 11 . 

Reference [5T] updates [33] and already includes 
KLOE's results. The recently released data of 
CMD-2 [23] and SND |23] are not yet included. 

The authors of [3j>] pioneered the idea of us- 
ing vector spectral functions derived from the 
study of hadronic r decays (see |37I38| for re- 
cent reviews) to improve the evaluation of the 
dispersive integral. Indeed, assuming isospin 
invariance to hold, the isovector part of the 
cross section for e + e~ — > hadrons can be cal- 
culated via the Conserved Vector Current rela- 
tions from r-decay spectra. The latest analy- 
sis with ALEPH |3J|, CLEO gU|, and OPAL gl] 
data yields a£ LO = 7110 (50) exp (8) rad (28) su{2) x 
10 _11 [3S|. Isospin-breaking corrections were ap- 
plied 02]- Information from r decays was also 

1 The evaluation of 1291 is not included as its result is being 
revised 1301 . 



included in one of the analyses of [31] , leading to 
a™° = 7027(47) e;cp (10) ra(i x l(T n . 

Although the precise CMD-2 e + e~ — > ir + ir~ 
data 22 are consistent with the corresponding r 
ones for energies below ~ 0.85 GeV, they are sig- 
nificantly lower for larger energies. KLOE's it + it~ 
spectral function confirms this discrepancy with 
the t data; on the contrary, the recent SND re- 
sults are compatible with them [21]. This dis- 
crepancy could be caused by inconsistencies in 
the e + e~ or t data, or in the isospin-breaking 
corrections which must be applied to the latter. 
Indeed, the mentioned disagreements between the 
e + e~ data sets need careful consideration. On 
the other hand, in spite of the agreement of the 
r data sets |38| , the question remains whether all 
possible isospin-breaking effects have been prop- 
erly taken into account [31143) . 

The hadronic higher-order contribution can be 
divided into two parts: a" HO = a" HO (vp) + 
a™°{\hl). The first term is the 0(a 3 ) contri- 
bution of diagrams containing hadronic vacuum 
polarization insertions [H|. Its latest value is 
< HO (vp) = -97.9(0.9) e;cp (0.3) rad x 10" 11 jSH], 
obtained using e + e~ annihilation data; it changes 
by ~ —3 x 10 -11 if hadronic r-decay data are used 
instead g5]- The second term, also of 0(a 3 ), 
is the hadronic light- by- light contribution; as it 
cannot be determined from data, its evaluation 
relies on specific models. In 2001 the authors 
of |46| uncovered a sign error in earlier evalua- 
tions of its dominating pion-pole part. Their esti- 
mate, based also on previous results for the quark 
and charged-pions loop parts gZ], is a^ HO (\bl) = 
80 (40) x 10~ n HE]. A higher value was obtained 
in 2003 including short-distance QCD constraints: 
a^ ff °(lbl) = 136(25) x 10" 11 g^. Further inde- 
pendent calculations would provide an important 
check of this result for a* HO (ibl), a contribution 
whose uncertainty may become the ultimate lim- 
itation of the SM prediction of the muon g— 2. 

4. Standard Model vs. Measurement 

The first column of Table 1 shows a^ M = a® ED + 
a* w +a% L °+a% H °. The values employed for a% LO 
are indicated by the reference on the left f |34| 
quotes two values, see Sec. 3); all a^ M values were 
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us 


116 591845 (69) 


235 (91) 


2.6 (3.0) 


E2 


116 591859 (90) 


221 (108) 


2.1 (2.5) 


EU 


116 591835 (69) 


245 (91) 


2.7 (3.1) 


El 


116 591855 (55) 


225 (81) 


2.8 (3.2) 


E3 


116 592 018 (63) 


62 (87) 


0.7 (1.3) 


G2 


116 591938 (54) 


142 (81) 


1.8 (2.3) 



Table 1: Standard Model vs. measurement. 



derived with a™°(M) = 136 (25) xl0~ n 49^. Er- 
rors were added in quadrature. The present world 
average experimental value for the muon g— 2 is 
a EXP = 116 592 080(60) x 10~ n (0.5 ppm) 0. 
The differences A = af, XF — af. M are listed in 
the second column of Table 1, while the num- 
bers of "standard deviations" (a) appear in the 
third one. Higher discrepancies, shown in angle 
brackets, are obtained if a" HO (lbl) = 80(40) x 
10- 11 HE] is used instead of 136 (25) x 10" 11 02]. 

5. Conclusions 

The discrepancies between recent SM predic- 
tions of a M and the current experimental value 
vary in a very wide range, from 0.7 to 3.2 a, 
according to the values chosen for the hadronic 
contributions. In particular, the leading-order 
hadronic contribution depends on which of the 
two sets of data, e + e~ collisions or r decays, are 
employed. The puzzling discrepancy between the 
7r + 7r~ spectral functions from e + e~ and isospin- 
breaking-corrected t data could be caused by in- 
consistencies in the e + e~ or r data, or in the 
isospin-breaking corrections applied to the latter. 
Indeed, disagreements occur between e + e~ data 
sets, requiring further detailed investigations. On 
the other hand, t data sets are in agreement, but 
their connection with the leading hadronic con- 
tribution to a p is less direct, and one wonders 
whether all possible isospin-breaking effects have 
been properly taken into account. Using e + e~ 
data, the SM prediction of the muon g— 2 deviates 
from the present experimental value by 2-3 a. 

The impressive results of the E821 experiment 
are still limited by statistical errors. A new ex- 
periment, E969, has been approved (but not yet 
funded) at Brookhaven in 2004 |4I50| . Its goal 
is to reduce the present experimental uncertainty 
by a factor of 2.5 to about 0.2 ppm. A letter of 
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intent for an even more precise g—2 experiment 
was submitted to J-PARC with the proposal to 
reach a precision below 0.1 ppm While the 
QED and EW contributions appear to be ready to 
rival these precisions, much effort will be needed 
to reduce the hadronic uncertainty by a factor of 
two. This effort is challenging but possible, and 
certainly well motivated by the excellent oppor- 
tunity the muon g~2 is providing us to unveil (or 
constrain) "new physics" effects. 
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